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Abstract. We study 'purely radiative' (i.e. divE = divH = 0) and geodesic perfect 
fluids with non-constant pressure and show that the Bianchi class A perfect fluids can 
be uniquely characterized — modulo the class of purely electric and (pseudo-)sphcrically 
symmetric universes — as those models for which the magnetic and electric part of 
the Weyl tensor and the shear are simultaneously diagonalizable. For the case of 
constant pressure the same conclusion holds provided one also assumes that the fluid 
is irrotational. 
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1. Introduction 

In a previous paper [24] we discussed 'purely radiative' (i.e. divE = divH = 0) 
irrotational dust spacetimes for which the magnetic part of the Weyl tensor is diagonal 
in the shear-electric eigenframe. We show that this analysis can be generalized to 'purely 
radiative' and geodesic perfect fluids. These are solutions of the Einstein field equations 

Rab ~ T^RQab = (p + P)u a u b + pg ab (p + p^0), (1) 

where u a is the normalized 4-velocity, p and p the density and pressure of the fluid and 
g ab the spacetime metric, and for which we impose the additional restrictions (1) that 
the flow is geodesic: u a = u a - b u b = and (2) that the electric and magnetic parts E ab 
and H ab of the Weyl tensor have vanishing spatial divergence. The spatial tensors E ab 
and H ab are hereby defined by 

E a b = C acb( iu c u d = E <ab> H ab = -s ac dC cd be u e = H <ab> , (2) 

where e abc = r\ abcd u d is the spatial projection of the spacetime permutation tensor rj abc d 
and S' <a 6 > = h a c h b d S( C( i) — \Scdh cd h ab is the projected, symmetric and trace-free part 
of S a b, with h ab = g ab + u a u b the spatial projector into the comoving rest space. The 
covariant spatial divergence and curl for tensors are defined by [14J: 

(divS) a = D b S ab , cm\S ab = e cd(a D c S b) d , (3) 

where the spatial part D a of the covariant derivative V a is given by 

D a A b ... = h a c h b d ---V c A d .... 
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Gravitational radiation is covariantly described by the nonlocal fields E ab (the tidal part 
of the curvature), which generalizes the Newtonian tidal tensor, and by H a t,, which has 
no Newtonian analogue [7j. As such, H a b may be considered as the true gravitational 
wave tensor, since there is no gravitational radiation in Newtonian theory. However at 
least in the linear regime, as in electromagnetic theory, gravity waves are characterized 
by H a i, and E^, where both are divergence-free but neither is curl-free [9]. For this 
reason spacetimes in which div-E = = divH have been termed purely radiative |17j . 

Although, when considering the gravitational dynamics of matter in the full non- 
linear regime, these conditions give rise to two independent chains of severe integrability 
conditions, it was remarked in [17] that both are satisfied by the Bianchi class A models. 
This motivated us to investigate purely radiative perfect fluids in their full generality. 

As a first step we have shown that in fact all the non-rotating dust models belong 
to Bianchi class A, provided we assume that they are not purely electric (in which case 
the term 'purely radiative' would be inappropriate anyway). In the purely electric case 
the exceptions were given by the Petrov type D spatially homogeneous and locally 
rotationally symmetric (LRS) metrics which are pseudo-spherically symmetric (and 
hence of Bianchi type ///), or which are spherically symmetric and belong to the 
Kantowski-Sachs family of dust models and hence do not admit a 3 dimensional isometry 
group acting simply transitively on the hypersurfaces of homogeneity. For technical 
reasons we also assumed that [a, H] = 0. Together with the fact that [a, E] = 0, which 
is a consequence of divH = 0, this implies that a a b, E ab and H ab can be simultaneously 
diagonalized (it is not necessary to assume that [E, H] = 0, as a degenerate shear 
eigenplane also must be a degenerate E eigenplane, see section 5). 

When, on the contrary, one considers fluids with non-constant pressure (but 
still with geodesic flow), the spatial gradient of p vanishes and the velocity becomes 
orthogonal to the p = constant hypersurfaces, such that the flow is irrotational [16| 118]. 
It is then a straightforward consequence of the assumptions divE = and [a, H] = 
that the spatial gradient of the matter density p vanishes too. As stated in [16] it 
follows that p, a and 6 are constants over the hypersurfaces of constant density and 
hence the resulting spacetimes satisfy the so called Postulate of Uniform Thermal 
Histories [H El El [20]. This seems to suggest that a local group of isometries 
exists, mapping the flow lines into each other (so called observational homogeneity). 
However not much progress has been made in this area and a detailed proof, or 
even a precise formulation of the conjecture that PUTH would lead to observational 
homogeneity, is still lacking. In fact it was suggested in [2] that any extra distinguishing 
mathematical property leading to observational homogeneity probably would have to be 
very complicated, involving at least third-order derivatives of the metric. An example of 
such an extra property is precisely the condition of being purely radiative, in the above 
sense of having divE = = divif, which indeed enables us to show that the resulting 
spacetimes are spatially homogeneous and, more particularly, of Bianchi class A when 
we exclude the 'non-radiative' purely electric subcases of section 15.31 
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2. Basic variables and equations 

We follow the notations and conventions of the Ellis-MacCallum orthonormal tetrad 
formalism [15]. Herein, the normalized 4- velocity u a plays the role of the timelike basis 
vectorfield of the tetrad {eo,e Q }. Greek indices take the values 1, 2 and 3 and refer to 
tetrad components with respect to e a . 

The basic variables in the formalism are 24 independent linear combinations of the 
Ricci rotation coefficients or, equivalently, of the commutator coefficients ^ d associated 
with the tetrad; these are the objects n a/3 and a a defined by = e^ e rf a + 8$ap — 5aas, 
the components Q of the angular velocity of the triad e Q with respect to the local 
'inertial compass' and the components it a , u a and 9 a p of the spatial kinematic quantities 
(acceleration vector, vorticity vector and expansion tensor respectively). We also use 
the decomposition Q a p = o a p + 9/35a/3, with a a p the trace- free shear tensor and 6 the 
expansion scalar. It turns out that there is a slight advantage in using instead of the 
variables n a! 3 and a a the connection coefficients n a , q a and r a related to the former by 

n a+1 „_! = (r Q + q a )/2, a a = (r a - q a )/2, n aa = n a+1 + n a _ x 

(these expressions have to be read modulo 3, so for example a = 3 gives n 12 = q?, + r 3 ). 
For sake of simplicity we define Z a , m a and S a as the spatial gradients of 8, p and 
a 2 = a ab a ab respectively: Z a = d a 9, m a = d a p and S a = d a a 2 . 

Together with the Einstein field equations, the basic equations are the Bianchi 
equations and the Ricci equations (part of which is equivalent with the Jacobi equations), 
which we have included in the appendix. 

3. Generalities 

From the Bianchi identity ( 1361) for geodesic perfect fluids (see appendix) and the 
assumption divH = we get, 

(div#) Q = DPHvp = ^ [E, a] = 0, 

which guarantees the existence of a frame wherein E a p and a a p are both diagonal and 
which will be used henceforth. 

Because all off-diagonal elements of E a p vanish we derive from ( 1321) that: 

^i(o"22 - 033) =02(033-011) = 3 (cr u - cr 22 ) = 0. (4) 

As a consequence the angular velocity Q a can be assumed to be 0: either all components 
automatically vanish, or, when e.g. the (e 2 , e^) plane is a shear eigenplane, then 
2 = 3 = and we can choose an extra rotation in this plane to make O x = 0. 
The latter operation leaves the frame then fixed up to rotations in the (e 2 , es)-plane for 
which the rotation angle (p satisfies do<p = 0. In fact, under a rotation about an angle <p 
in the (e 2 , e3)-plane the quantities O transform as follows: 

1; 2 , 3 — > Q\ — dop>, 2 cos ip — O3 sin ip, 2 sin + O3 cosy?. 
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The Raychaudhuri equation (1431) and the conservation law (1341) give us the evolution 
of respectively the expansion and the density, while (|3T|) gives the evolution of the 
components a aa : 

3<Vii = (0-22 + °33 ~ 2cr n ~ 2o- n 6») - 3E U , 

3d a 22 = (o-g 3 + o\ x - 2o 22 - 2o 22 0) - 3E 22 , (5) 

39 cr 3 3 = {<A.\ + °22 - 2(T 33 ~ 2o 33 6») ~ 3£ 3 3- 

The spatial derivatives of obtained by the (Oa)-Einstein field equations (T44l) 

and by the off-diagonal components of ( l28l) : 

01011 = +ri(cr u - cr 22 ) - gi(cr n - 0-33), 

5 2 crii = - § ^2 + 92(0-22 - era )- #13, (6) 
<9 3 crn = -~Z 3 + r 3 ((7 U - cr 33 ) + i7 i2 

and similar expressions by cyclic permutation of the indices. Note that the divE = 
condition relates the off-diagonal parts of H a p to m a by (1351) : 

mi = 3H 2 3(a 33 - <t 22 ), 

m 2 = 3^13(0-11-0-33), (7) 
m 3 = 3H 12 (a 22 - on). 

Using (1341) we can now act with the commutator [do, d a ] = — (o QQ + \9)d a on p to obtain 
the relation: 

4 

<%m a = -m a (a aa + -6 1 ) - (p + p)Z a , a = (1,2,3). (8) 

Propagating ([7]) along the fluid flow and substituting ([5]) and (jHJ) gives then the 
following expressions: 

3(o u - a 22 )d Q H 12 = (6o n - 6o 2 2 2 - 6E 22 - 3E 33 + 2o 22 # - 2a u 9)H 12 + (p + p)Z 3 
3(o 22 - a 33 )d H 23 = (6o 2 2 2 - 6o 3 2 3 - 6E 33 - 3E U + 2o 33 ^ - 2a 22 9)H 23 + (p + p)Z 1 (9) 
3(033 - OnH#i3 = (6o 3 2 3 - 60^ - 6E U - 3E 22 + 2a n - 2a 33 9)H 13 + (p + p)Z 2 . 

Combining the equations above with the spatial information in the Bianchi and 
Jacobi equations allows one to obtain explicit algebraic expressions for all the directional 
derivatives of the variables p, p,9,a aa , Z a , E aa , H a/ 3,r a ,q a and n a . Although the 
integrability conditions for the resulting Pfaffian system are polynomial and although 
the complicated nature of these restrictions seems to suggest that solutions can only 
be spatially homogeneous, the size of the polynomials involved and the large number of 
different subcases left to be investigated, make it impossible so far to draw any general 
conclusions. Therefore we will restrict in the present paper to one of the main subcases, 
namely [o, H] = 0. As mentioned in the introduction this assumption implies that H a p 
is diagonal in the (o, £')-eigenframe (in the case of degenerate shear (13~T1) implies that 
the eigenplanes of o a/3 and E a p coincide). From ([7]) and ([8]) or ([9]) one obtains then 
m a = and Z a = (regardless of the degeneracy of o Qj a), giving a first hint that the 
corresponding spacetimes might be spatially homogeneous indeed. To prove this we first 
construct some information about the remaining part of the magnetic tensor. 
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Using the [d , d a ] commutators on 9, we obtain from (}3"Uj) and Z a = m a = d a p = 
that S a = d a a l 3~ / <J l3 ' y = and hence, using §6§, 

n{cru ~ 022) 2 - gi(o"n - 0-33) 2 = 0, 

r 2 (cr 2 2 ~ cr 33 ) 2 - q 2 (a 22 ~ cr u ) 2 = 0, (10) 
r 3 (o"n - cr 33 ) 2 - q 3 ((T33 - CX22) 2 = 0. 

We also need the spatial derivatives of E aa , which we can deduce from ( 1351) and 
the off diagonal part of ( 1381 ) : 

d x E n = E 11 {n-2q 1 )-E 22 {r 1 + q 1 ), 

d 2 E u = (E 22 -E u )q 2 , (11) 
d 3 E u = (E u - E 33 )r 3 . 

We will treat the two cases of non-degenerate and degenerate shear now separately. 
Note that by ([5]) and fl33|) a a p 7^ 0, as otherwise the solutions would be FLRW. 



4. Non degenerate shear 

As the frame is invariantly defined, we first aim to proof that r a and q a are zero and 
that all spatial derivatives vanish. We first look at the evolution of a 2 : from (1301) it is 
easy to check that 

4 

d a 2 = — 9a 2 - 2a 3 + 2E-a, (12) 
3 

where a 3 = a ia affa a ^ and E • a = E a p<j afi . 

Hence, acting with the commutator [<9 , d a ] on a 2 we find: 

d S a = ^9S a - ^a 2 Z a - 2d a a 3 - 2d a (E ■ a) - S a a aa (13) 

and hence, as Z a = S a = 0, d a a 3 + d a (E ■ a) = 0. Working out the latter componentwise 
and subsituting (TTTT) and (Q, results in three algebraic expressions. For a = 1 we get 

ri(a 2 2-^n) [4^22+2^33-30-33 (a 2 2-o-ii)]-gi(o-33-o- 1 i)[2 J E 22 +4 J E33-3o-22(o- 3 3-crii)] = 0.(14) 

Eliminating r\ from ( fTUl) and ( fT4l) . we find 

91(0-33 - (Tn){<J-z2 - <?n)x = ( 15 ) 

with 

2 

X = gt-E 1 !^^ - 033) + ^((Tss - an) + £33(0-11 - a 22 )} - (a u - o 22 )(p 22 - a 33 ){p 33 - a u ) 

The whole calculation can be repeated for the cases a = 2, 3 so that we finally have (as 
the shear is non-degenerate) 

HX = Q2X = Q3X = 0. (16) 

If x 7^ this implies with ffTUj) that r a = q a = 0, while when x — one has dix = 
and hence, using ([6]) and ( JTTT) . 

^1(0-33 - 0-n) 3 + n(o- 2 2 - 0-n) 3 = 0. (17) 
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Eliminating qi from (TTTj) and (ITUla) and using tr(a Q , / 3) = again leads to r a = q a = 0. 

This means that a a = and n a p is diagonal. It is easy to check that 
all spatial derivatives of the rotation coefficients now vanish and hence [6] we 
obtain the non-degenerate spatially homogenous Bianchi class A models, i.e. types 
J, II, VI , VIIq, VIII and IX. 

5. Degenerate shear 

From the degeneracy of the shear it is easily verified that E a p must be degenerate too 
(I3T]) and we may assume 

(Tap = diag(-2cr, a, a) 
E aP = diag(-2£,£,£). 

For dust we can discard the case of vanishing E as a consequence of the 
non-existence of anti- Newtonian universes [25], but for non-constant pressure purely 
magnetic perfect fluids, which are necessarily of Petrov type I or D, do exist (we 
don't consider the conformally flat these uniquely reduce under the present 

assumptions to the FLRW models, see [ID]): for type I the only possibility is the 
Bianchi VIq metric discussed in [23], while for type D the spacetime is locally 
rotationally symmetric class III [23] . As the fluid is non- rotating the solutions are then 
given precisely by the Bianchi type VIII or IX Lozanovski-Aarons-Carminati (LAC) 
metrics [12j [13], [TT] , with a possible degeneracy to Bianchi type II (see below). 

Henceforth we take E non zero (both the type I and type D purely magnetic 
models above will reappear as E = special cases). Assuming that H a/3 is diagonal in 
the (a, E)-eigenframe, H af3 = diag(-H 22 - H33, H 22 , H 33 ), CGI) implies 

ri - qi = r 3 = q 2 = 0, (18) 

while the (23)-component of (1291) shows that r% + q\ = and hence r\ = q% — 0. 

At this point the only variables having possibly non-vanishing spatial derivatives 
are q 3 ,r 2 and ni, n 2 ,n s . With the simplifications obtained so far, the (ll)-Einstein 
equation ( H5i) reduces to 

= l -p + l -e{a - l -6) + ^a 2 -E- n 2 n 3 . (19) 

Eliminating now E a p from ( 137|) for a = (3 = 1 and a = (3 = 2, we obtain 

= #n (-2m + n 2 - 2n 3 ) + F 33 (-4m - n 2 - n 3 ), 
which in combination with ( J33l) leads to 

{ri! + n 2 + n 3 ){n 2 - n 3 ) = 0. (20) 

We will discuss the three cases that follow from fl20l separately; i.e. n 2 7^ n 3 , n 2 = n 3 7^ 
and n 2 = n 3 = 0. Notice that when n 2 = n 3 the necessary and sufficient conditions for 
the fluid to be locally rotationally symmetric are automatically satisfied [8]. 
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5.1. n 2 -n 3 ^0 

By ( |20l) we have n\ + n 2 + n 3 = 0. 

From (jHJ) we get the evolution of n a , namely 
d m = a(n 2 + n?) - \n x 6 
don 2 = -o(n 2 + 3n 3 ) - \n 2 6. 
Herewith we obtain, propagating rii + n 2 + n 3 along e , 

= d (ni + n 2 + ^3) = — 3(n 2 + n 3 )<7 
and thus rii = and n 2 = —n 3 7^ 0. Applying <9 3 to ffl9l) we find <9 3 n 2 = 0. On the other 
hand, calculating d 3 n 2 from the (12)-Einstein equation ( H6l) and the Jacobi equation 
(|39|) . results in d 3 n 2 = q 3 {n 3 — n 2 ). Hence we see that 

q 3 {n 3 - n 2 ) = 

and thus q 3 = 0, as n 3 7^ ri 2 by assumption. It follows now that also r 2 = as <9ig 3 = r 2 rz 2 
by ( l46l) and ( J39l) . So r a and g a all become zero and we have 

n af3 = diag(0, -n 2 ,n 2 ) (n 2 ^ 0) 

a Q = 0. 

This implies the vanishing of the spatial derivatives of all rotation coefficients (d a = 0) 
and we obtain a spatially homogeneous universe of Bianchi class A, type VIq. Examples 
are the metrics without rotational symmetry in [1]. 



5.2. n 2 = n 3 j^0 

From ( fl9i) and the (33)-Einstein equation ( )45l) we can calculate 

8 2 r 2 =p + 3a 2 - \d 2 -n\- 2n 2 m + q 2 3 + r\ + d 3 q 3 . (22) 
At this stage n a/3 and a Q are given by 





' 2n 2 q 3 /2 


r 2 /2 







n a p = 


q 3 /2 nx + 712 







r 2 /2 




r 2 /2 


ni + n 2 




-qs/2 



Now 9i?2 2 = by the Jacobi equation (|39l) . while 9 2 n 2 = and 9 3 n 2 = follow by 
combining respectively the (31)-Einstein equation (H6l) and the (12)-Einstein equation 
f|46|) with the Jacobi equation fl39|) . so the spatial derivatives of n 2 vanish. We now try 
to obtain the Bianchi A condition a a = 0: remembering that we still have a rotational 
degree of freedom left (namely rotations in the (23)-plane about an angle ip satisfying 
d ip = 0), we use the following transformation formulas for the quantities n a/ 3 and a a : 
n a /3 -> n r a/3 , a a -> a r a with 

n r u = n n , 

n i2 = ^[ c os((p)(d 2 ip + q 3 ) - sxn((p)(d 3 (p + r 2 )], 
n i3 = ^[ cos (v)(d 3 ip + r 2 ) + sin(v?)(<9 2 v3 + g 3 )] 
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and 

n 23 = a[ = a r 2 - n r 13 = a r 3 + n\ 2 = 0. (24) 
A rotation making a a = and n a p diagonal can then be obtained by d 2 ip — ~Q3 
and d^ip = —r 2 , the integrability condition of which leads to di^p = — ^ + n i + n 2, with 

K = p + 3cr 2 -\d 2 + n\. (25) 

One can check that the integrability conditions for the resulting system of pde's 
d cp = 0, 

c/i <p = h ni + no, 

^ 2n 2 2 ' (26) 

<9 3 ^ = - r2 

are identically satisfied. For the rotated variable n\ we find ri\ + n 2 = all spatial 
derivatives now vanish and we obtain the spatially homogeneous LRS Bianchi class A 
spacetimes of types II (K = 0), VIII (K < 0) and IX (K > 0). An integrable subcase 
arises when p = 3n 2 — 8 (a — 8/3) — Qa 2 : the electric part of the Weyl tensor is then 
and we obtain the purely magnetic LAC metrics [T2l [T3| [TT] , for which the Bianchi 
type is VIII or IX, reducing to type II for the p = p/5 Collins-Stewart metric [3]. As 
Ea/3 = the label 'purely radiative' is however not fully appropriate for these models. 

Of course the above reasoning breaks down when n x = n 2 = 0, a case which will be 
dealt with below. 



5. 3. n 2 =n 3 = 

Within the previous class the n 2 = n 3 = models are exceptional in the sense that 
they are 'purely electric' (by fl33l) one has H a/3 = 0) and again cannot be termed 'purely 
radiative'. For completeness however we present the full details of the resulting special 
cases. 

As in (5.2) one shows that a rotation exists under which r 2 , ni and d^q^ become 
zero. By (|T9l) this implies g| + K — 0, where K is the same quantity as defined by (1251) . 
The rotation is now determined by a solution of the following system: 

d ip = 0, 

ft*> = (27) 

d 2 ip = cos((p)y/-K - q 3 , 

d 3 (p = - sm((p)V-K - r 2 

for which again the integrability conditions are identically satisfied. 
When K < 0, one obtains 








V-K/2 


" 







n a p = 


V-K/2 


























-V-K/2 
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which is a LRS spatially homogenous Bianchi class B type III model [B|. 

When K = the rotation fl27j) results in a a = 0, n a( g = 0, which gives a LRS 
spatially homogenous Bianchi class A type / model (equivalent to Bianchi class A type 
VII [6J). 

A problem arises when K > 0: qs, n\ and ri2 can still be made zero by a rotation, 
but it is no longer possible to choose a frame in which all spatial derivatives identically 
vanish. However the Cartan equations can now be integrated easily and one finds the 
Kantowski-Sachs perfect fluids, with metric given by 

ds 2 = -dt 2 + Q(t) 2 dr 2 + P(t) 2 [d6 2 + sin 2 (#)# 2 ] 

and P(QP - PQ) + P(QP - PQ) + Q = 0. As is well known these metrics do not 
admit a 3-dimensional isometry group acting simply transitively on the t = constant 
hypersurfaces. 

6. Conclusion 

We have shown that the Bianchi class A perfect fluid models with non-constant pressure 
can be uniquely characterized as geodesic perfect fluid spacetimes which are purely 
radiative in the sense that the gravitational field satisfies divif = divE = 0, under 
the assumption that also the magnetic part of the Weyl tensor H a/ 3 is diagonal in the 
shear-electric eigenframe (i.e. [H, a] = 0). The only possible exception arises in the 
purely electric Petrov type D case of degenerate shear, where the allowed solutions are 
the Kantowski-Sachs perfect fluids or the pseudo-spherically symmetric Bianchi type III 
models (K = p + 3cr 2 — |# 2 > resp. K < 0). The same conclusion holds for constant 
pressure provided one also assumes that the fluid is irrotational. 

The next logical step is to investigate the case where H, E and a are not 
simultaneously diagonalizable: it is hoped that here new classes of inhomogeneous and 
cosmologically interesting solutions might turn up, or, alternatively, that one would be 
able to demonstrate the remarkable result that the Bianchi class A spacetimes are the 
unique (modulo the LRS exceptions above) 'purely radiative' ones. Aside from the fact 
that such a characterization of the Bianchi A models would be quite neat, we feel that 
a further investigation of restrictions on divi? and divif and of the role played by these 
third order quantities in general relativistic perfect fluids might shed a new light on the 
issues of PUTH and observational homogeneity. 
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7. Appendix 

A first set of basic equations consists of Ricci and Bianchi equations [7J, which we 
will use in a more formal sense and therefore give in their covariant form (assuming 
u a = u a = 0): 

Ricci equations 

curiae -H ab = 0, (28) 
2 

o- <a b> + -9o- ab + a c<a o b> c + E ab = 0, (29) 

9 + \e 2 + <7 a6( 7 a6 + X -{p + 3p) = 0. (30) 

In their orthonormal tetrad form (with E a p and o a p diagonal) the equations ([281 
and (|29|) read respectively 

#n = -<Vn + g(o-| 2 + 4, - 2(7^ - 20o-u), (31) 
fii(*22 - 033) = (32) 



and 



#11 = (^22 - 0"ll)^3 + (0"33 - 0"ll)«2, 

1 1 1 (33) 

#12 = 2^3(0-11 - 022) + -(cr 33 - crn)r 3 + -(0-33 - 0-22)^3 

and similar equations obtained by cyclic permutation of the indices. 

Bianchi equations 

P = -(p+p)0, (34) 

div£ a = [a, H] a + ^D a p, (35) 

dwH a = -[a,E} a , (36) 

E <a b> - cm\H ab = -9E ab + 3a c<a E b> c - - (p + p)a ab , (37) 

H <a b> — cwc\E ab = —9H ab + 3a c<a H b> c , (38) 
where the spatial dual of the commutator of tensors, [S, T] a is defined as [S, T] a = 

ZabcO J-d ■ 

A second set of equations consists of those Jacobi and Einstein equations which we 
use explicitly in our calculations and which we present in their orthonormal tetrad form 
(using the eigenframe of a ab and E ab ), but now with the extra simplification u a = u a = 
and Q a = (see section [3]). Most of these equations appear in triplets which can be 
obtained from each other by cyclic permutation of the indices. For each equation we 
give only one representant. 
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Jacobi equations 



(J2- 


J'4): 


(J'8- 


jlO) : 


(Jll 


-Jl3) 


(Jl4 


-jl6) 



<9i(n 2 + n 3 ) + cV 3 + <9 3 g 2 = (39) 
(ri - qi)(n 2 + n 3 ) + r 3 r 2 - q 3 q 2 , 

d (n - qi ) - d 1 (a 11 - h) = -(n - qi )(<Jn + ^), (40) 

d {n 2 + n 3 ) = (2a u - -6){n 2 + n 3 ), (41) 

d (r 3 + q 3 ) + d 3 (a u - or 22 ) = -(a 33 + ^6>)(r 3 + g 3 ). (42) 



Einstein equations 



(EinOO) : <9 £ = -[(<t u + ^) 2 + (<7 22 + ^) 2 + (<7 33 + ^) 2 ] - ^(p + 3p), (43) 

2 1 3 

(EinOa) : -drf - d&u = -(n + q 1 ){a 22 - cr 33 ) - -(n - gi)<7n, (44) 

(Einaa) : - <9 ((Jii + -0) - di(n - q x ) + d 2 q 2 - d 3 r 3 = ^(p - p) (45) 

+ 92(^2 - q 2 ) - r 3 (r 3 - q 3 ) - r\ - q\ + 2n 2 n 3 + 6>(o- n + ^0), 
(Eina/3) : - <9ir 2 + <9 2 <?i - d 3 (ni - n 2 ) = qi(r 2 + q 2 ) + r 2 (n + 2^) (46) 

- (r 3 - gaXwi - n 2 ) + 2n 3 (r 3 + g 3 ). 

Commutator relations 

Finally we have the commutator relations, which can be written succinctly as 
[di, d 2 \ = q 2 di + r\d 2 + (m + n 2 )<9 3 . 
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